Multiband tight-binding theory of disordered A^B^C semiconductor quantum dots: 
Application to the optical properties of alloyed Cd^Zn^^Se nanocrystals 
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Zero-dimensional nanocrystals, as obtained by chemical synthesis, offer a broad range of appli- 
cations, as their spectrum and thus their excitation gap can be tailored by variation of their size. 
Additionally, nanocrystals of the type A^Bi-^C can be realized by alloying of two pure compound 
semiconductor materials AC and BC, which allows for a continuous tuning of their absorption and 
emission spectrum with the concentration x. We use the single-particle energies and wave functions 
calculated from a multiband sp 3 empirical tight-binding model in combination with the configuration 
interaction scheme to calculate the optical properties of Cd^Zni-zSe nanocrystals with a spherical 
shape. In contrast to common mean-field approaches like the virtual crystal approximation (VCA), 
we treat the disorder on a microscopic level by taking into account a finite number of realizations 
for each size and concentration. We then compare the results for the optical properties with recent 
experimental data and calculate the optical bowing coefficient for further sizes. 

PACS numbers: 73.22.Dj, 78,67.Hc, 71.35.Cc, 71.23.-k 
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I. INTRODUCTION 

As a special realization of semiconductor quantum dots 
(QDs), photochemically stable semiconductor nanocrys- 
tals (NCs) can act as an alternative to organic molecules 
in a broad range of applications, such as light-emitting 
devices (LEDs)ji^ laser applications^ and especially for 
biological fluorescence labelingj^£ 

One common approach to directly tune the pho- 
toluminscence (PL) or absorption spectrum of these 
nanocrystals in the ultraviolet-visible range is by varia- 
tion of their size. This approach has extensively been 
used with common II- VI and III-V compound semi- 
conductor NCs and these systems (mostly using CdSc, 
ZnSe, CdS and core-shell-combinations of those) are 
part of a large variety of experimental^ 11 as well as 
theoretical^— studies. 

In order to obtain emission at shorter wavelengths (500 
down to 400 nm) for common CdSe NCs, QDs with a 
very small diameter of less than 2 nm have to be pro- 
duced^ which corresponds to less than 4 conventional 
lattice constants a. In this size regime, the surface pas- 
sivation and size control is very difficult and one is left 
with a low PL efficiency* 2 ^, When using materials with 
larger bulk band gaps like ZnSe or CdS, the diameter of 
the NCs must be enlarged and the efficiency is then low- 
ered by the relatively weak confinement. An alternative 
is the usage of alloyed NCs of the type A x Bi_ x C with 
either varying or fixed size, as the emission wavelength 
can additionally be tuned by variation of the concentra- 
tion x. Such alloyed NCs have for example been realized 
for Cd^Zni^S^ and Cd^Zni.^Se^^ Particularly, the 
latter Cd^Zni-^Se NC system allows for the coverage 
of the whole visible spectrum with a fixed diameter of 
approximately 3 nm. 

Similar to the behaviour of most mixed bulk semicon- 
ductors, the energetic position of the single particle gap 



(given by the difference e\ — hi of the lowest bound elec- 
tron state and the highest hole state) as well as of the 
absorption/emission lines of Cd^Zni-^Se NCs shows a 
pronounced bowing behaviour as a function of the con- 
centration x. The deviation from linear behaviour is com- 
monly described using a single bowing parameter b. The 
concentration dependent energetic position of the single- 
particle gap, respectively the spectral line of the alloyed 
NCs is then given by 

E(x) = xE AC + {l-x)E BC -x(l-x)b, (1) 

where -Bac and -Ebc are the corresponding quantities of 
the pure AC and BC material. In the case of binary 
bulk alloys, also a non-parabolic dependance of the band 
gap on x is observed when the pure materials AC and 
BC have considerably different lattice constants^ which 
should directly carry over to the QD properties. Never- 
theless, a fit of experimental data by the ansatz (JTJ is still 
common and mostly possible. Even in the case of a per- 
fectly parabolic behaviour, the bulk band gap _Eg Ulk (a;), 
the single-particle gap E^ c (x) and the position of the 
PL peak E Fh (x) will of course not show the same bow- 
ing behaviour. 

While the single-particle states and energies them- 
selves are not observed directly in optical measurements, 
their proper calculation is of crucial importance from a 
theoretical point of view. They give the foundation for a 
proper many-particle approach for the optical properties, 
like the Hartree-Fock method^ 7 - and the configuration in- 
teraction (CI) schemei 27 ' 28 The CI uses a finite subspace 
spanned by the single-particle wave functions of the QD 
system to construct the many-particle states. 

For the calculation of the single-particle spectrum of 
unalloyed zero-dimensional QD systems, a broad spec- 
trum of computational methods with different levels of 
sophistication can be used, like effective mass approxima- 
tions^iir— multiband k • p-modelsj22r— empirical pseu- 
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dopotential models (EPM)i£iii2& M and empirical tight- 
binding models (ETBM) ^^^ii 

Within the ETBM approach, one can either start from 
atomic orbitals localized at the sites of the crystal struc- 
ture (LCAO)p^ or directly use Wannier-like orbitals lo- 
calized at the sites of the Bravais lattice, i. e. on each 
unit cell. This defines the so-called effective bond-orbital 
model (EBOM)4£ri£ As each unit cell in a semiconductor 
compound of the type A x Bi_ x C is either occupied by the 
AC- or BC-Material, the spatial resolution of the EBOM 
is especially suitable for these mixed systems. The basic 
idea in ETBM/EBOM approaches is to deduce a set of 
equations for the TB matrix elements in terms of a finite 
number of bulk parameters (taken from experiments or 
first-principle calculations) in order to adequately repro- 
duce the bulk band structure. The same matrix elements 
then serve as input in the calculations for the nanostruc- 
tures. 

In order to find a suitable approximation to calcu- 
late the single-particle states for alloyed systems within 
a TB model, it is possible to map the microscopically 
disordered system on an effective system without disor- 
der by either using a concentration-dependent interpola- 
tion of relevant parameters (virtual crystal approxima- 
tion, VCA) or more sophisticated Green function meth- 
ods (coherent potential approximation, CPA). The VCA 
is known to vastly underestimate even the bulk bowing 
in most cases4^ While the CPA can in fact give good 
results for bulk crystals j 50 ' 51 it becomes much more com- 
plex and less reliable when applied to low-dimensional 
systemsi 52 ' 53 

An alternative for bulk as well as QD systems is to 
explicitely perform calculations for an ensemble of mi- 
croscopically distinct realisations, which is computation- 
ally very costful, but exact in the disorder and can re- 
produce bowin g 49 ' 54 and line broadening effects^ The 
single-particle states can be used in CI calculations to 
determine the many-body states and energies and, there- 
fore, measurable quantities like the absorption/emission 
spectrum. 

As it has been shown recently^, the treatment of the 
EBOM with disorder on a finite supercell yields excellent 
results for bulk Cd^Zni-^Se. In order to test the validity 
of this approach when applied to quantum dots, we use 
our model to calculate the concentration-dependent emis- 
sion spectrum of violet- to orange-emitting Cd^Zni-^Se 
alloy nanocrystals. Such systems have experimentally 
been realized by means of a cation exchange reaction by 
Zhong et al. in Ref. 24, and have been optically charac- 
terized by their absorption and PL spectra in the same 
work. 

This paper is organized as follows. In SecHH we 
present the theoretical approach for our calculation of 
the electronic and optical properties of the Cd x Zni_ K Se 
alloy nanocrystals. It is based on an empirical tight- 
binding model with disorder on a microscopic level, and 
the calculation of dipole and Coulomb matrix elements 
from the single-particle states. Section IIII Al contains a 



TABLE I: Material parameters for zincblende CdSe and ZnSe, 
as used in Ref. IlO . 



Parameter 


Description 


CdSe 


ZnSe 


a (A) 


lattice constant 


6.078 


5.668 


E g (eV) 


band gap 


1.76" 


2.82 


AE vh (eV) 


valence band offset 


0.22 


0.00 


Xi (eV) 


X-point energy of CB 


2.94 


4.41 


XI (eV) 


X-point energy of HH/LH band 


-1.98 


-2.08 


XI (eV) 


X-point energy of SO band 


-4.28 


-5.03 


A so (eV) 


spin-orbit splitting 


0.41 


0.43 


m c (m ) 


CB effective mass 


0.120 


0.147 


7i 




3.33 


2.45 


72 


Kohn-Luttinger parameters 


1.11 


0.61 


73 




1.45 


1.11 



"taken from Ref.l56l. 



detailed discussion of the resulting single-particle spec- 
trum, while the following section IIII Bl discusses the in- 
fluence of the disorder on the dipole and Coulomb ma- 
trix elements. In Sec. IIII Cl we compare the resulting 
concentration-dependent optical spectrum to experimen- 
tal results for a fixed nanocrystal size, while Sec. IIII Dl 
deals with the size-dependence of the bowing of the opti- 
cal lines and the question whether the bowing parameter 
can be approximated by the bulk limit for larger sizes. 



II. THEORY 

A. Effective bond-orbital model for pure bulk 
semiconductors 

The optical properties of direct band gap semiconduc- 
tors stem from an s-like conduction band (CB) and three 
p-like valence bands, the heavy hole (HH), light hole (LH) 
and split-off (SO) band. Consequently, we use one s- and 
three p-orbitals per spin direction, localized on the sites 
R of the fee Bravais lattice underlying the zincblende 
crystal structure: 

|Ra) , ae{s \,p x t,Py t,Pz t,sl,p x l,P v UPz 1} ■ 

(2) 

The TB matrix elements of the bulk Hamiltonian 
H hulk are then given by-- 

E^' = (Ra\H hulk \R,'a') . (3) 

In order to fit the bulk band structure to the bulk band 
gap iJ g , the conduction band effective mass m e , the spin- 
orbit splitting A so , the conventional lattice constant a, 
and additionally to the X-point energies of the conduc- 
tion band Xi c , the HH and LH bands X$ v (which are 
degenerate at the Brillouin zone boundary), and the 
SO band X^ v , we use the parametrization scheme as 
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FIG. 1: (Color online) Realization scheme and geometry of 
CdxZni-icSe nanocrystals by cation exchange reaction, as 
realized by Zhong et al. (following Ref. 24) . The resulting 
nanocrystals have a diameter of slightly above 3 nm. 



presented by Loehr in Ref. |48l The necessary num- 
ber of free parameters is achieved by the restriction of 
the non-vanishing hopping matrix elements up to sec- 
ond nearest neighbors and the usage of the two-center- 
approximation^ 5 - The choice of material parameters for 
CdSe and ZnSe is presented in Tab.|TJ This set is, with 
one exception, identical to the one used in Ref. |4j] and 
has been proven to be reliable in similar TB calculations 
for pure CdSe and ZnSe quantum dots. It originally 
stems from Refs. l58ll59l augmented by X-point energies 
from RefM 



B. Application to ensemble of alloyed nanocrystals 

The CdzZni^Se nanocrystals from Ref. [24] are spher- 
ical in shape. Their diameter before the alloying by a 
cation exchange reaction is d = (3.2 ± 0.3) nm for the 
CdSe NCs and d = (3.1 ± 0.3) nm for the initial ZnSe 
NCs, as determined by transmission electron microscopy. 
Their geometry and realization is schematically shown in 
Fig.[T] Due to the spherical shape and their flexible sur- 
rounding, NCs without core-shell-structure show almost 
no strain^i 

In analogy to previous EBOM calculations for unal- 
loyed zero-dimensional nanostructureS) 57 ' 62 ! 63 we model 
a spherical NC with the diameter d, where the localized 
orbital basis of Eq. © is located on each site of the un- 
derlying fee lattice. As each of the sites belonging to the 
Cd x Zni_ x Se NC will either be occupied by the diatomic 
CdSe or ZnSe basis, the probability for the respective 
species is given by the concentrations x and 1 — x when 
assuming uncorrelated substitutional disorder. Accord- 
ingly, we will use the TB matrix elements, Eq. ([3]), of the 
pure CdSe or ZnSe material for these lattice sites. For 
hopping matrix elements between unit cells of different 
material, we use the arithmetic average. 

The confinement potential at the NC boundaries is 
assumed to be infinite and thus incorporated by set- 
ting the hopping to zero at the surface. This common 
approac h 18 ' 19 i 41 i 64 corresponds to a perfect surface passi- 
vation by hydrogen or ligand molecules. Previous studies 
have shown that the effect of a finite barrier of appropri- 
ate height on the optical properties is negligible in the size 
regime which is considered in the present paper ^ The 



valence band offset (VBO) between CdSe and ZnSe is 
incorporated by shifting the respective site-diagonal ma- 
trix elements by the value A£ v b = 0.22 eV4^ A desired 
finite number of single-particle electron and hole states 
\t/jf ), \ipi) and energies ef, e^ 1 can finally be computed by 
numerical diagonalization of the resulting Hamilton ma- 
trix. 

It should be pointed out that within the EBOM 
scheme, each lattice site of the A^Bi-^C alloy can be 
unambigously assigned to the site-diagonal TB matrix el- 
ements and the corresponding material parameters for ei- 
ther AC or BC. In microscopic ETBM calculations which 
resolve the anion and cation sites, disorder on an atomic 
scale prevents this, as each anion of the type C will lo- 
cally be surrounded by a different number of A or B 
cationsi 49 i 65 This makes an assignment of the diagonal 
elements of the anions to the band structure of either 
AC or BC impossible and will lead to an effectively more 
coarse-grained resolution as in the case of the EBOM. 

Previous EBOM supercell calculations on bulk 
Cd x Zni_ x Se^ have shown that the bowing behaviour 
safely converges when the finite ensemble contains at 
least 50 microscopically different configurations per con- 
centration. The same number has been found to be suffi- 
ciently reliant in the present calculations, i. e. the numer- 
ical results for the average single- particle and PL energies 
can be reproduced with the same degree of accuracy (10 
meV) as the input parameters from Tab.Q] In addition to 
the configurational disorder, we also include the effects 
of concentrational disorder, 55 - by not enforcing the overall 
concentration per configuration to equal the local prob- 
abilities for finding CdSe or ZnSe. This is crucial for the 
correct reproduction and prediction of experimental find- 
ings, because the experimentally determined AC or BC 
concentrations refer to the proportions in the ensemble 
and not in the single configuration. 



C. Coulomb and Dipole Matrix Elements, 
Many-body Hamiltonian 

The many-body Hamiltonian in terms of creation and 
annihilation operators e\(h\) and e^/ij for electrons 
(holes) in the single-particle state \i/)f) with en- 

ergy ef(e^) reads 

H = H +H C + H D . (4) 
It consists of the diagonal single-particle part 



the Coulomb interaction part 



He = 



t t 

ijkt e i e j e kei 



ijkl 



(6) 
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and a part which describes the coupling to an external 
field E in dipole approximation: 

H »=Y, (<V>i|eoE • r|$) e^- + h.c.) , (7) 



where en is the electron charge. The present form of 
He omits the contributions from electron-hole exchange 
terms, as they are very small compared to the other ma- 
trix elements (see Ref.[(36| for further details). According 
to the discussion in Ref. HH, the Coulomb matrix ele- 
ments Vijki can be approximated by 

v$ = £ E ( c r«) * (46) * 4&<& v& r') , 

(8) 



RR' act' 



where the asterisk denotes the complex conjugation, 



V(R - R) 



47re £ r |R- R'| 



for R ^ R' 



(9) 



where p = E/|E| is the light polarization vector. Strictly 
speaking, the above given expressions r and dfj give only 
the envelope part and neglect the spatial dependence in- 
side the unit cells, as the latter will only give minor con- 
tributions to the overall value. We refer to Ref.|63| for a 
more detailed discussion on this topic. 

Even in unalloyed QD systems, band mixing effects 
prevent the use of selection rules based on the total an- 
gular momentum, as it is not a good quantum number 
any more.— Additionally, the configurational disorder in 
our alloyed system will in general break any spatial sym- 
metry of the NC Hamiltonian, so that each realization 
will show unique dipole and Coulomb matrix elements, 
which underly a certain distribution for the ensemble. In 
contrast to mean-field simulations, which will preserve 
the symmetry of the potential, this can be accounted for 
in our finite-ensemble approach by explicitly calculating 
the single-particle spectrum, dipole and Coulomb matrix 
elements for each concentration and configuration. 



and 



uc J uc 



d 3 rd 3 r'- 



47T£o£ r |r — r' 



Vb • (10) 



Here, V uc = a 3 /4 is the volume of the unit cell of the 
fee lattice. According to Ref.H3, the calculation of the 
on-site integral V(0) can be done quasi- analytically by 
expansion of the Coulomb interaction in terms of spher- 
ical harmonics. The coefficients c RtJ! are numerically ob- 
tained as entries of the i th eigenvector of the TB matrix 
and are related to the corresponding one-particle wave 
function by 



Ra 



(r) , 



(11) 



where </>R a (r) denotes the a- type effective orbital local- 
ized at the Bravais lattice site R, see Eq. (0. Note that 
for finite sized nanocrystals in solution, the proper choice 
of the dielectric constant e r is more complicated than in 
the case of overgrown self- assembled QD structures. We 
will adress this topic explicitly in the next section. 

In order to obtain the selection rules and oscillator 
strengths, one has to calculate the matrix elements d^j 1 = 
e (V , fl r l 1 /'^) of the dipole operator e r using the EBOM 
wave functions ipi (r) . In the dipole Hamiltonian, Eq. (0 , 
we make use of the approximatio n 19 ' 68 ' 69 



E |Ra)R(Ra| 



(12) 



Ra 



which represents an approximation for the position op- 
erator r in a TB approach. Then, the TB dipole matrix 
elements dfj = p • dfj 1 explicitly read 



RR' 



( c rL) c R ' l Q ,p-R<5 RR '(W, (13) 



D. Treatment of charge screening 

Due to the small diameter of the NC system under con- 
sideration, finite size effects have to be taken into con- 
sideration when calulating the Coulomb matrix elements 
with an appropriate, site-averaged dielectric constant e r . 
It has been shown before that £ r will be significantly 
smaller than its bulk value for small QDs (d < 5 mn), 
due to a combination of finite size effects and the energy 
gap increase as induced by the confinement i 1 5 ' 7 1 1 72 Be- 
cause the Coulomb matrix elements scale directly with 
l/e r , a careful choice of the dielectric response is crucial 
for a proper reproduction and prediction of experimental 
spectra, especially as the Coulomb interaction increases 
significantly for smaller NCs. 

In the present paper, we will use a similar approach 
as Lee et al. in Ref. [l8|. The general idea is to use a 
dielectric function which is dependent on both the diam- 
eter d of the quantum dot and the separation r between 
two particles. The result is a combination of the mod- 
ified Thomas-Fermi model proposed by Resta^ 3 - for the 
seperation dependence and a generalization of the Penn 
model 7 -^ for the size dependence. It has first been used in 
this combination in Ref. [13 for unalloyed NCs, where this 
approach is discussed in detail. Explicitly, this dielectric 
function reads 

e S) C ( d ) Vol [ sinh (<?( r - r o)) +qr], r <r 



ErM)= le»), 



where 



(jjbulk + A )2 

[£ g NC (d) + A 



el c (d) = l + ( £ £*-l) 



2 ' 



r > r , 
(14) 



(15) 



Here q = (4/7r) 1 / 2 (37r 2 n ) 1 / 3 is the Thomas-Fermi wave 
vector, containing the valence electron density uq = 
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32/a 3 of the zincblende structure. In our TB approach, 
the particles are represented by respective portions of the 
charge density as obtained by the of Eq. (JTTJ) and 
thus r re |R — R'|. This is on the level of a monopole 
approximation and consistent with the afore-mentioned 
treatment of the Coulomb and dipole matrix elements, 
ro is a critical radius, which has to be determined by the 
solution of sinh (qro)/qr = e^(d) (see Ref.ll7l). e^ lk is 
the high-frequency dielectric constant of the bulk mate- 
rial, and A is basically the bulk exciton binding energy, 
which can be obtained from A = E 2 — -Eg Ulk < 0, where 
Ei is the energy of the first pronounced peak in the bulk 
absorption spectrum. Following experimental data ) 56 i 75 
we will use A(CdSe) = 15 meV and A(ZnSe) = 20 meV 
in our calculations. 

Especially for the zincblende modification of CdSe, 
the values for e^ lk differ in the literature to a cer- 
tain degree, partially due to the fact that cubic CdSe is 
metastable and the experimental growth of high-quality 
samples is very difficult and often depends on the spe- 
cific substrate^ When using standard literature val- 
ues, (e.g. Refs. 56.75), this can lead to the spurious 
choice of £^ lk (CdSe)/e^ lk (ZnSe) w 1, which is not suit- 
able for calculating the properties of the present mixed 
Cd x Zni_ x Se systems, as the dielectric response should 
become larger when the bulk band gap becomes smaller. 
In order to obtain a consistent parameter set, we use a 
version of the empirical Moss model for II- VI semicon- 
ductors, 7 £ where 



_bulk 



lk ■ 



(16) 



In this approximation, k is an empirical parameter, 
with k « 108 eV for II-VI-semiconductors. In combi- 
nation with the parameter set from Tab. HI this yields 
e^ lk (CdSe) « 7.83 and e^ lk (ZnSe) sa 6.19. The latter is 
in good agreement with established experimental values 
for ZnSeM 

At this point, we already anticipate that the use of the 
resulting single-particle gaps E^ c (d) for different diam- 
eters d of the pure or mixed NCs in Eq. (IT5|) always gives 
a screening radius r$ in the range of 1.5-2.5 A. Thus, it 
does not exceed the nearest-neighbour distance of the un- 
derlying fee lattice and according to Eq. (TH1) we will use 
e r = e^ c (d) in the Coulomb matrix elements V(R — R') 
with R ^ R'. For the site-diagonal Coulomb interaction 
V(0), we follow the guidelines of Ref.[l|| In this work, 
Lee et al. performed Monte Carlo calculations with TB 
orbitals on similar, unalloyed Si and CdSe NCs, with the 
result that the effectice screening of the on-site Coulomb 
integrals is approximately half the long-range screening 
given by £^p(d). In large QD systems, the resulting value 
of the Coulomb matrix elements is not sensitive to the 
exact value of V(0), due to the long-ranged nature of 
the Coulomb interaction. We carefully checked that this 
statement still remains valid in the finite systems under 
consideration, which is also in accordance with the results 
of Ref.QJ. 
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FIG. 2: (Color online) Resulting scaling factor l/e^ c (d) from 
the dielectric function for pure CdSe and ZnSe nanocrystals as 
function of the diameter d. For the sake of comparison, the 
corresponding scaling factor l/e^ lk for the CdSe and ZnSe 
bulk materials is also given. 



For the mixed Cd^Zni-^Se NCs, the differentiation be- 
tween CdSe and ZnSe sites for the size-dependant dielec- 
tric function £ r (r, d) would not serve any purpose. Be- 
cause any truly microscopic dielectric response function 
would additionally be dependent on the positions of the 
inducing and responding charges, the commonly used e T 
is an effective, site-averaged quantity anyway. The scal- 
ing factors l/e^(d) for pure CdSe and ZnSe nanocrys- 
tals as function of the diameter d are given in Fig.[5J 
together with the corresponding bulk values l/ej^ lk . As 
the NC scaling factor for the pure materials only differs 
slightly, especially for the experimentally realized size 
(d w 3 nm) of Ref.iU, we can safely use a concentra- 
tionally averaged dielectric constant, 



r (r, d,x) 



xe 



CdSe 



(r,d) + (l-x)ef nSc (r,d), (17) 



when calculating the Coulomb matrix elements of the 
Cd^Zni-^Se nanocrystals. 



III. RESULTS FOR MIXED NANOCRYSTALS 

A. Single-particle energies 

In this section, we will discuss the single-particle spec- 
trum of the Cd x Zni_ x Se nanocrystals as obtained from 
our TB calculations on a finite ensemble. As already 
discussed in detail in Section [II Bl we will assume spher- 
ical quantum dots with a diameter d between (3.1 ± 0.3) 
nm (pure ZnSe) and (3.2 ± 0.3) nm (pure CdSe). The 
discretization on the fee lattice allows for a spatial res- 
olution of the confinement potential of half the conven- 
tional lattice constant, so that the experimental size can 
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FIG. 3: Distribution of the single-particle energies of the six lowest electron levels ei_6 and six highest hole levels hi-e of the 
CcUZni-^Se nanocrystals for TV = 50 realizations. The levels of gray give the relative frequency of the respective values. For 
the disordered NCs (0 < x < 1), the levels show a strong broadening, especially towards the Zn-richer concentrations. Each 
of the states e, and hi shows an additional twofold Kramers degeneracy due to time reversal symmetry. All energies are given 
with respect to the valence band edge of ZnSe. 



TABLE II: Single-particle energies of the six lowest electron 
levels ei_6 and six highest hole levels hi-e for the unalloyed 
CdSe and ZnSe nanocrystals. Each of the states shows an 
additional twofold Kramers degeneracy. All energies are given 
with respect to the valence band edge of ZnSe. 



CdSe NCs 




ZnSe NCs 




ei (eV) 


2.377 


ei (eV) 


3.269 


e 2 (eV) 


2.703 


e 2 (eV) 


3.657 


e 3 (eV) 


2.707 


e 3 (eV) 


3.659 


e 4 (eV) 


2.707 


e 4 (eV) 


3.659 


es (eV) 


3.014 


es (eV) 


4.045 


e 6 (eV) 


3.014 


e 6 (eV) 


4.045 


hi (eV) 


0.089 


hi (eV) 


-0.114 


hi (eV) 


0.089 


hi (eV) 


-0.114 


hz (eV) 


0.082 


h z (eV) 


-0.119 


h 4 (eV) 


0.082 


hi (eV) 


-0.119 


h 5 (eV) 


0.026 


h 5 (eV) 


-0.139 


h 6 (eV) 


0.026 


h 6 (eV) 


-0.139 



be modelled within the error boundaries. We account for 
50 different microscopic realizations per concentration x, 
including both the effects of configurational and concen- 
trational disorder. 

As the single-particle energies and symmetry proper- 
ties of the corresponding eigenfunctions of unalloyed II- 
VI nanocrystals have been extensively discussed in the 
literaturef ^ 15 ' 17 ' 20 ' 21 i 36 ' 41 we will here concentrate on 
the influence of the disorder on a finite number of eigen- 
states and energies. 

Figure [3] shows the distribution of the single-particle 



energies of the six lowest electron levels ei_6 and six high- 
est hole levels hi-e per spin direction of the Cd^Zni-^Se 
NCs for N — 50 realizations, where the levels of gray 
give the relative frequency of the eigenvalues. Each of 
these levels is twofold degenerate due to the Kramers 
time-reversal symmetry of the system. For the pure QD 
systems (x = 1 and x = 0), the lines show additional 
degeneracies that stem from the spatial symmetries of 
the underlying potential. The corresponding eigenener- 
gies are additionally given in Tab.HIl and show an ex- 
cellent agreement with the results of Ref.[4l|, where the 
spectrum and its degeneracies of pure CdSe nanocrys- 
tals have been discussed by means of a fully microscopic 
Scp^ empirical tight-binding model and have furthermore 
been compared to experimental results by Alperson et 
al& Figure [3] also clearly shows that the energy levels of 
the mixed systems (0 < x < 1) are considerably broad- 
ened due to the microscopic disorder. The broadening is 
more prominent for lower concentrations x of the small 
band gap material (CdSe). Especially in case of the hole 
levels for x < 0.5, the disorder- induced variation of the 
eigenenergies exceeds the initial level splitting, so that 
an unambiguous assignment to the initial level order is 
not possible anymore. This alone renders the informative 
value of the results of simple mean-field approaches like 
the virtual crystal approximation questionable. 

We additionally want to examine whether the virtual 
crystal approximation, where the TB matrix elements are 
obtained as a concentration-dependent linear interpola- 
tion, can at least reproduce the behaviour of ensemble- 
averaged quantities and choose the single-particle energy 
gap Eg C = ei — hi as an example. The results are 
depicted in Fig. 2) The ensemble-averaged energy gap 
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FIG. 4: (Color online) Average energy gap Eg C,av for each 
concentration x of the CdrcZni-^Se nanocrystals, obtained 
from the values depicted in Fig. [3] The black line gives the 
ensemble-averaged results for TV = 50 realizations per con- 
centration, and the red circles the VCA results. Additionally, 
the energy gap e\{x) — h\(x) of each microscopically distinct 
configuration is given by the black triangles. The VCA gives 
an eronneous linear behaviour over the whole concentration 
range. 



£ g NC ' av (i) = l/NJ2?=i (4( x ) ~ h i( x )) shows a distinct 
bowing behaviour as a function of the concentration x. 
Clearly, the results from the VCA approach cannot re- 
produce any bowing behaviour and give a perfectly linear 
dependence of the gap on the concentration. Further- 
more, even the influence of small fluctuations from the 
pure cases (e. g. x — 0.1 and x — 0.9) is vastly underesti- 
mated. For the sake of comparison, we have also added 
E % E (x) = e\(x) — h\(x) for each microscopically distinct 
configuration, which corresponds to the energy gaps of 
the single nanocrystals in the ensemble. 

In order to concentrate on the comparison with optical 
measurements, we will not further analyze the bowing be- 
haviour of the single particle properties here, but proceed 
to calculate the excitonic spectra from the results of the 
TB calculations on the finite ensemble with microscopic 
disorder. 





FIG. 5: (Color online) Ensemble-averaged modulus square 
of the dipole matrix element d% h for two concentrations x of 
the Cd^Zni-jjSe nanocrystals for TV = 50 realizations (each 
normalized with respect to the maximum value). Degener- 
ate states are merged into a common field. In the pure case 
(x = 0), clear-cut selection rules can be obtained. Due to the 
loss of the spatial symmetries in the alloyed case (x — 0.5), 
all degeneracies (besides the essential Kramers doublets) are 
lifted and the strict selection rules are relaxed. 
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FIG. 6: (Color online) Distribution of the Coulomb matrix 
element V1331 for each concentration x of the Cd^Zni-^Se 
nanocrystals for TV = 50 realizations. In the pure cases (x = 
and x = 1), this matrix element gives the main contribution 
to the redshift of the lowest excitonic absorption line. In the 
alloyed cases (0 < x < 1), the values show a broad distribu- 
tion. 



B. Dipole and Coulomb matrix elements 

From the single-particle wave functions, the dipole ma- 
trix elements and Coulomb matrix elements V££, can 
be calculated in the framework given in the sections III CI 
and III Dl for each concentration and configuration. 

The ensemble-averaged modulus square of the dipole 
matrix element d%r for two concentrations x of the 
Cd x Zni_ x Se nanocrystals, each normalized with respect 
to the maximum value, is given in Fig. [5] Here, d% h is 
the projection on the (arbitrarily chosen) [100] direction. 



In the pure case, (x = 1 and x = 0), clear-cut selection 
rules can be obtained. The energetically lowest allowed 
transition for the x — 1 and x — (not shown here) 
case is e\ — h^i, in agreement with previous results from 
theory and experiment ! 24 ! 77 ' 78 Due to the loss of the spa- 
tial symmetries in the alloyed case (exemplarily shown 
for x — 0.5), all degeneracies, the essential Kramers dou- 
blets left aside, are lifted and the strict selection rules are 
relaxed. 

In order to study the influence of the disorder on the 
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Coulomb matrix element, we consequently choose the el- 
ement Vi3 31 , which will give the main contribution to 
the redshift of the lowest excitonic line in the sense of a 
diagonal approximation (see Ref. 70 for further details) 
in the pure case. The distribution of the values of Vx33i 
for each microscopic configuration is depicted in Fig. [6] 
for each concentration. In the pure cases, the relative 
frequency is of course given by a ^-distribution, as there 
is only one possible configuration. In the alloyed cases 
(0 < x < 1), the values show a more or less broad distri- 
bution. 



C. Calculation of optical spectrum and comparison 
to experiment 

In order to calculate the excitonic absorption/emission 
spectrum for each concentration x of the Cd^Zni-^Se 
nanocrystals, we start from the calculated single-particle 
wave functions, dipole and Coulomb matrix elements for 
each microscopic configuration and perform CI calcula- 
tions in order to obtain the many-particle states. Since 
we are only interested in the energetically lowest transi- 
tions, we include four electron and four hole states per 
spin direction. Then, the absorption/emission lines can 
be calculated by usage of Fermi's golden rule, 



: S{E l -E f ±hu), (18) 



where \^>i) (|\&/)) is the initial (final) many-particle state 
with energy E$ (£7), as obtained by diagonalization 
of the Hamiltonian (j3J). For more details, we refer to 
Refs.lH and [lO. The resulting spectrum for the NC en- 
semble with concentration x is then obtained by the su- 
perposition of all N = 50 spectra. 

Although there is indication that the Stokes shift in II- 
VI nanocrystals is either induced by the electron-hole ex- 
change interaction ; 77 ' 79 or exciton-acoustic phonon cou- 
pling^ there also exists experimental evidence that the 
underlying mechanism is more complex*^ As our approx- 
imation for the many-body Hamiltonian does not contain 
corresponding terms anyway, we will not differentiate be- 
tween absorption and emission lines in our calculations 
from here on. Instead, we slightly modified the TB pa- 
rameters in order to reproduce the energetically lowest 
emission line in the pure cases (x = 1 and x = 0) as 
obtained from the experiment, which ensures the correct 
boundary values for the investigation of the bowing be- 
haviour. This is the usual approach in similar TB calcu- 
lations for bulk systems ; 49 ' 65 and here necessary to elim- 
inate further influences of e. g. the finite size resolution of 
both theory and experimental characterization, ambigu- 
ities in the input parameter set (especially for CdSe, as 
discussed in Sec. Ill Dl) and the temperature dependence 
of the optical spectra. In retrospect, the bowing be- 
haviour turned out to be insensitive to this procedure 
(not shown). 
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FIG. 7: (Color online) Optical spectrum for x — 0.5 for the 
CcUZni-^Se nanocrystals as obtained by the ETBM calcula- 
tions on the finite ensemble. The two quasi-broadened peaks 
each develop from distinct, multiply degenerate emission lines 
in the pure case, due to the disorder in the single-particle en- 
ergies, dipole and Coulomb matrix elements. 



In Fig. we examplarily show the optical spectrum 
for x — 0.5. The clearly visible two "broadened" peaks 
each develop from two distinct, multiply degenerate emis- 
sion lines in the pure case when the disorder is intro- 
duced. Ultimately, this peak structure is the consequence 
of the disorder in the single-particle energies, dipole and 
Coulomb matrix elements as depicted in the Figs. [3J [3] 
and IH which introduces a different shift and line height 
for each configuration. From the envelope of these broad- 
ened peaks we can read off the energetic position of the 
lowest transition of the finite ensemble in dependence of 
x. The 50 realizations per concentration turn out to be 
sufficient to rcproducibly determine this position with an 
accuracy of 10 meV (which corresponds to the resolution 
of the input parameters, see Tab.QJ for the NC size under 
consideration. 

The concentration-dependent results are summarized 
in Fig. [51 where the lowest excitonic transition of the 
Cd K Zni_ :r Se NC ensemble is given as a function of the 
concentration x. Additionally, we give the experimen- 
tal results from the optical characterization by Zhong et 
al. for the system under consideration. Obviously, there is 
good agreement with our theoretical curve, albeit all ex- 
perimental values are slightly shifted to higher energies. 
We calculate a bowing of b = 0.66 ± 0.08 eV, while the 
experimental results give the smaller value b — 0.5 ±0.1 
cV. The relatively large error shows that the common 
assumption of a parabolic dependence on the concentra- 
tion is only of limited applicability. Nevertheless, both 
values agree within the error boundaries, which is espe- 
cially noteworthy when comparing to bulk Cd^Zni-^Se, 
where a broad range of bowing values between 6 = 
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FIG. 8: (Color online) Energy of the lowest excitonic transi- 
tion of the ensemble as a function of the concentration x, 
as obtained from the ensemble-averaged spectrum for the 
d « 3.2 nm Cd^Zni-^Se NCs. The black squares give the 
tight-binding ensemble results , while the red triangles give 
the experimental results from Zhong et al, as given in Ref. 
[H. We calculate a bowing of b = (0.66 ±0.08) eV, while their 
results give b = (0.5 ± 0.1) eV. 



and b = 1.26 eV is reported. 65 ' 82 " — Furthermore, a com- 
parison with similar TB-calculations and corresponding 
experimental results for unalloyed nanocyrstals (e. g. as 
m Refs-EHl) reveals the remarkably good concordance 
of our results with the experiment. 



D. Calculations for further sizes and reproduction 
of bulk limit 

In this section, we finally want to discuss the de- 
pendence of the bowing parameter b on the nanocrys- 
tal diameter d. As we found a good agreement for the 
Cd ^Zni—^Se nanocrystals realized by Zhong et al. with 
d ~ 3.2 nm, we calculated the optical properties for larger 
sizes, where the TB formalism should be more accurate. 
Furthermore, we simulated the "bulk" limit d — » oo by 
using a sufficiently large supercell (4000 lattice sites) with 
periodic boundary conditions (see Ref.Hi] for more de- 
tails). Again, N = 50 realizations were used for each 
concentration. The results are given in Fig. [9] and show 
that the bowing parameter reasonably approaches the 
bulk limit value when the diameter is increased. 

Obviously, the influence of the finite size of the nanos- 
tructure on the nonlinear behaviour of the optical gap 
is rapidly decreasing. This raises the question, whether 
the bowing behaviour for larger sizes could properly be 
incorporated by a modified VCA approach, as suggested 
by di Carlo in Ref.l85l. Such an approach has been used 
in several calculations for mixed QD systems, like e. g. 
in Ref.[H for the calculation of excitonic properties of 



FIG. 9: (Color online) Size dependence of the bowing pa- 
rameter b. Here, d is the concentrationally-averaged diameter 
of the Cd^Zni-zSe NCs. As before, we used N = 50 re- 
alizations. The values approaches a bulk limit, which has 
been calculated with a sufficiently large supercell with peri- 
odic boundary conditions. 




FIG. 10: (Color online) Average energy gap £^ c,av for 
each concentration x of the Cd^Zni-^Se nanocrystals for the 
largest calculated size of d « 6.8 nm. The black squares give 
the ensemble-averaged results for TV = 50 realizations, the 
red circles the VCA results and the green triangles the results 
from the modified VCA, where the bowing parameter of the 
bulk limit d — > oo has been incorporated. The modified VCA 
results coincide well with the ensemble-averaged gap for the 
Cd-richer concentrations. 



small embedded lens-shaped In^Gai-^N/GaN QDs with 
an eight-band k • p-model. 

Here, the diagonal TB matrix elements are weighted 
nonlinearly with the bulk bowing parameter, such that 
the correct bowing is empirically reproduced in the bulk 
limit. The results for the average single-particle gap 
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^■NC.av are depicted in Fig. [10] In contrast to the lin- 
ear VCA approach, the modified VCA can reproduce 
the average behaviour of the single-particle gap over 
the concentration x quite well, especially for large Cd- 
concentrations. Nevertheless, if one is interested in the 
influence of small fluctuations of the small band gap ma- 
terial, the modified VCA still does not give accurate re- 
sults, as the influence of the finite size on the energy 
levels and the disorder obviously spoils the transferabil- 
ity of the bulk bowing to the QD energies. For smaller 
diameters, the deviation between the ensemble-averaged 
and the modified VCA results becomes more prominent 
(not shown). 

The advantage of the modified VCA is the signifi- 
cantly lesser effort in the numerical calculations once the 
proper bulk bowing is accessible (we already discussed 
the broad range of the experimental bowing values, be- 
cause of which we calculated the bulk bowing ourselves 
within the same TB approach). In the modified VCA 
case, only one calculation per concentration is necessary, 
so that it can surely be a feasible method if one is only 
interested in the coiicentrationally averaged properties of 
large systems. But although there were efforts to connect 
the bulk bowing parameter b to the TB parameters for 
the pure systems within the VCA, they fail to satisfac- 
torily reproduce the bowing in accordance with experi- 
mental results without free parameters. 87 ' 88 So the main 
drawback of this approach remains the fact that the mod- 
ified VCA does not calculate the bowing itself but only 
reproduces it from an external experimental or calculated 
bulk bowing, in contrast to our TB calculations with mi- 
croscopic disorder. Furthermore, no broadening effects 
can be simulated in the modified VCA. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we presented a theoretical approach for 
the determination of the electronic and optical properties 
of quantum dots consisting of binary compound semi- 
conductor alloys of the type A x Bi_a;C and applied it to 
spherical Cd^Zni-^Se nanocrystals. We used a multi- 
band empirical tight-binding approach with a basis set 
localized on the sites of the underlying Bravais lattice, 
which allows for an exact treatment of substitutional dis- 
order on the microscopic scale. We discussed the result- 
ing single-particle energies of a finite ensemble of alloyed 
nanocrystals as a function of the concentration x. Fur- 
thermore, we showed by comparison to results obtained 



by the virtual crystal approximation that such simple 
mean-field approaches cannot adequately simulate the in- 
fluence of the alloying and the disorder. From the single- 
particle energies and wave functions, we calculated the 
excitonic spectrum of a finite ensemble of 50 distinct re- 
alizations per concentration by means of a configuration 
interaction scheme. Special attention was paid to a care- 
ful choice of the dielectric constant, in order to ensure a 
consistent treatment for the calculation of the Coulomb 
matrix elements of the mixed system. 

The nonlinear concentration dependence of the spec- 
tral lines was then compared to experimental results as 
obtained by Zhong et al. for NCs with a diameter of 
about 3 nm, yielding a very good agreement between ex- 
periment and our theory. 

Finally, we presented results for larger nanocrystals 
and found that the bowing behaviour reasonably con- 
verges to the bulk limit. By using the bowing parameter 
b of the bulk limit as an input parameter, one can perform 
a modified VCA calculation. We showed that even the 
modified VCA is still of only limited applicability even in 
the case of fairly large (« 7 nm diameter) quantum dots. 

In contrast to similar approaches for alloyed systems 
like Ref. l86l which treat the disorder within a mean- 
field framework like the virtual crystal approximation, 
our simple TB supercell approach is able to satisfacto- 
rily reproduce the experimental findings for the alloyed 
system without additional free parameters. In case that 
the computational resources are available, our approach 
can easily be applied to further alloyed QD systems, ei- 
ther realized by chemical synthesis or epitaxy. For sys- 
tems where lattice strain plays an important role, e. g. 
embedded QDs with a high lattice mismatch, a proper 
calculation of a strain field might have to be additionally 
included for each realization, which will of course increase 
the calculational efforts. 
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